Prank's constant in the hexatic phase 
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Using video-microscopy data of a two-dimensional colloidal system the bond-order correlation 
function Ge is calculated and used to determine the temperature-dependence of both the orienta- 
tional correlation length ^ in the isotropic liquid phase and the Frank constant Fa in the hexatic 
phase. Fa takes the value 72 /n at the hexatic <-> isotropic liquid phase transition and diverges at 
the hexatic <-> crystal transition as predicted by the KTHNY-theory. This is a quantitative test of 
the mechanism of breaking the orientational symmetry by disclination unbinding. 

PACS numbers: 64.70.Dv,68.35.Rh,82.70.Dd 



The theory of melting in two dimensions (2d) de- 
veloped by Kosterlitz, Thouless, Halperin, Nelson and 
Young (KTHNY-theory) suggests a two-stage melting 
from the crystalline phase to the isotropic liquid. The 
first transition at temperature T m is driven by the disso- 
ciation of thermally activated dislocation pairs into iso- 
lated dislocations breaking the translational symmetry 
0, 0] ■ The fluid phase directly above T m still exhibits 
orientational symmetry and is called the hexatic phase. 
It may be viewed as an anisotropic fluid with a six-fold 
director [3J, |4| which is characterized by a finite value 
of Frank's constant Fa , the elastic modulus quantifying 
the orientational stiffness. At the second transition at 
Ti > T m , the dissociation of some of the dislocations 
into free disclinations destroys the orientational symme- 
try. Now, the fluid shows ordinary short-range rotational 
and positional order as it is characteristic of an isotropic 
liquid. 

Following an argument given in 0, , T m and Ti can 
be estimated using the defect interaction Hamiltonian 
Hd between a pair of disclinations (d — disc) and a pair 
of dislocations (d = disl) which for both defect pairs 
and at large distances goes like Hd ~ Qjlogr with the 
dimensionless strength parameter Cd depending on the 
defect type. Defect dissociation is completed at a tem- 
perature where the thermally averaged pair distance (n) 
diverges. Evaluating this expression for Hd one gener- 
ally finds divergence if Cd = 4. The unbinding condition 
Cd = 4 translates into lim T ^ T - j3JC(T)a^ — IQtt for dislo- 
cation pairs ((3 = l/fc^T, eto is lattice spacing) and into 
lim T ^ T - (3Fa(T) = 72 /it for disclination pairs, where 
/C is the Young's modulus of the crystal. Connecting 
thus the defect pair unbinding condition to the two tran- 
sition temperatures Ti and T m , two expressions are ob- 
tained that summarize the microscopic explanation of the 
KTHNY theory for two-stage melting. 

In this Letter we study the temperature-dependence of 
Frank's constant of a 2D system in the hexatic phase. We 
first determine the hexatic — > isotropic fluid transition 
temperature T and then check if Frank's constant takes 
the value 72/tt at T it thus testing the KTHNY theory 



and its prediction that disclination unbinding occurs at 
Ti. In addition, we analyze the divergence behavior of 
the orientational correlation length at T and of Frank's 
constant at T m . 

Different theoretical approaches invoking grain bound- 
ary induced melting 0, [|| or condensation of geometrical 
defects 0] suggest one first order transition. How- 
ever some simulations for Lennard-Jones systems indi- 
cate the hexatic phase to be metastable 0|. The 
transition in hard-core systems seem to be first-order 
probably due to finite-size effects Simulations with 
long-range dipole-dipole interaction clearly show second 
order behavior Experimental evidence for the hex- 
atic phase has been demonstrated for colloidal systems 
111 111 IH 111 El, in block copolymer films 0,1111, 
as well as for magnetic bubble a rray s and macroscopic 
granular or atomic systems (H Hi E El . Still 
the order of the transitions is seen to be inconsistent. 
The observation of a phase equilibrium isotropic/hexatic 
[ItI |2~0| and hexatic /crystalline 0] indicates two first 
order transitions. In our system we find two continuous 
transitions. 

The experimental setup is essentially the same as in 
[27|. Spherical and super-paramagnetic colloids (diam- 
eter d = 4.5 pm) are confined by gravity to a water/air 
interface formed by a water drop suspended by surface 
tension in a top sealed cylindrical hole of a glass plate. 
The field of view has a size of 835 x 620 pm 2 containing 
typically up to 3 • 10 3 particles (out of 3 * 10 5 of the whole 
sample). A magnetic field H is applied perpendicular 
to the air/water interface inducing in each particle a 
magnetic moment M = ■ This leads to a repulsive 
dipole-dipole pair-interaction with the dimensionless in- 
teraction strength given by V — f3{p Q / Att){xH) 2 (tt p) 3 / 2 . 
Here \ is tb_ e susceptibility per colloid while p is the 
2d particle density and the average particle distance is 
a = Xj^fp. The interaction strength can be externally 
controlled by means of the magnetic field H; it can be 
interpreted as an inverse temperature and is the only 
parameter controlling the phase behavior of the system. 
For each T the coordinates of the colloids are recorded 
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FIG. 1: Structure factor S(q) of our colloidal system at three different inverse temperatures F corresponding to the isotropic 
liquid (r = 52.4), the hexatic phase (r = 59.6) and the crystalline (r = 61.0) phase. The central cross is an artifact of the 
Fourier-transformation. 



via video-microscopy (resolution of particle position 
dr = 100 nm) and digital image processing over a period 
of 1 — 2 h using a frame rate of 250 ms. 

To set the stage we first visualize in Fig. {Q the three 
phases and their symmetries by plotting the structure 
factor 



S(q) 
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as calculated from the positional data of the colloids for 
three different temperatures. Here, a, a' runs over all 
N particles in the field of view while () denote the time 
average over 700 configurations. In the liquid phase, con- 
centric rings appear having radii that can be connected 
to typical inter-particle distances. The hexatic phase, on 
the other hand, is characterized by six segments of a ring 
which arise due to the quasi long-range orientational or- 
der of the six- fold director j2^. In the crystalline phase 
the Bragg peaks of a hexagonal crystal show up with a 
finite width that is due to the quasi long-range character 
of the translational order. 

To quantify the six-fold orientational symmetry the 
bond-order correlation function 



G 6 (r) = (V(r)V*(0)) 



(2) 



is calculated with ijj{r) — J2j e 6t9i ^ r \ Here the sum 
runs over the Nj next neighbors of the particle i at po- 
sition f and 9ij (f) is the angle between a fixed reference 
axis and the bond of the particle i and its neighbor j. 
() here denotes not only the ensemble average which is 
taken over all N(N—l)/2 particle-pair distances for each 
configuration (resolution dr = 100 nm) but also the time 
average over 70 statistically independent configurations. 



KTHNY theory predicts that 
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isotropic: short range , 



776 < 1/4 and takes the value 1/4 right at T = T t . All 
three regimes can be easily distinguished in Fig. show- 
ing Ge(r) for a few representative temperatures. Note, 
that Gq(0) is not normalized to 1. 

We next fit Ge(r) to r^ 6 and e~ r l^ & to extract r\% 
and £ 6 . The fits are performed for radii r/a £ {0..20} 
29]. To check for the characteristics of the orientational 
correlation function, the ratio of the reduced chi-square 




FIG. 2: Orientational correlation function Ge(r) as func- 
tion of the inverse temperature T in a log-log plot. From top 
to bottom: Three curves for the crystalline phase showing 
the long-range orientational order (Zim r _,ooG6(r) 7^ 0), two 
curves showing the quasi long-range order of the hexatic phase 
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and three curves showing the short-range 



order typical of the isotropic liquid (Ge(r) ~ e 
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FIG. 3: Quantitative test for the long distance behavior of 
Ge,(r). For Xalg/Xexp < 1 the algebraic decay fits better. 



the characteristic of the divergence. At T = 58.0 there is 
a huge finite size effect indicating that £e is much larger 
than the field of view. However, inside the hexatic phase, 
^6 is no longer well defined as the decay is algebraic. We 
fit our data to eq. © in the range 49 < T < 57.5 and find 
the critical exponent v = 0.5 ± 0.03 and Tj = 58.9 ± 1.1, 
a value which due to the finite-size effect is slightly larger 
than Ti obtained from Fig. Q. 

The exponent tjq is related to Frank's constant Fa 



% (T) 



18fc s T 



(4) 



So the critical exponent rj 6 (Ti) — 1/4 corresponds to 
(3FA{Ti) = 72 /tt at the hexatic «-» liquid transition. This 
quantity is plotted in Fig. (@Jd). Indeed, Fa crosses the 
value 72/7T at Ti = 57.5 ±0.5 exactly at that temperature 
which in Fig. has been independently determined to 
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exponential (xtxp) n ^ is shown in Fig. J2J as a function 
of r for three different measurements. For melting, a 
crystal free of dislocations was grown at high T and then 
r was reduced in small steps. For each temperature step 
the system was equilibrated 1/2 h before data acquisition 
started. This was done at different densities: melt A with 
average particle distance of a = 11.8 pan and meltJl with 
a = 14.8 fim containing 3200 respectively 2000 particles 
in the field of view. The measurement denoted freeze 
in Fig. © (a = 11.8 pro) started in the isotropic liquid 
phase and T was increased with an equilibration time of 
1 h between the steps. For Xai g /Xex P > 1 an exponential 
decay fits better than the algebraic and vice versa for 
Xalg/Xexp < 1- We observe in Fig. © that the change 
in the characteristic appears at Ti = 57.5 ± 0.5. This 
value is the temperature of the hexatic <-> isotropic liquid 
transition. 

In the vicinity of the phase transition, approaching 
Ti from the isotropic liquid the oricntational correlation 
length ^6 should diverge as 0, 



£ 6 (r)~exp (- 



i/r - l/r, 



U7) 



(3) 



with b a constant and v = 1/2. This behavior is observed 
in Fig. £ 6 indeed increases dramatically near Ti — 

57.5 ± 0.5 irrespective of whether the system is heated or 
cooled. Before discussing this feature we first address the 
finite size effect. To this end, we have computed Ge(r) 
and ^6 for subsystems of different size ranging from 720 x 
515 /ito 2 , 615 x 405 //to 2 , 505 x 300 /xto 2 , 400 x 190 pm 2 
to 390 x 80 yjrv? . The resulting data-points are plotted as 
triangles in Fig. (@J and belong to the black filled squares 
which they converge to. No finite size effect is found for 
r < 56, but a considerable one at T = 56.9 close to Ti 
where we obviously need the full field of view to capture 
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FIG. 4: Correlation length (a) and Frank constant Fa (b) 
as a function of the inverse temperature. ^6 diverges at Fi and 
Fa at F m . In between the system shows hexatic symmetry. 
The solid lines are fits to eq. (J^J and JSj resulting in critical 
exponents v = 0.5 ± 0.03 and v = 0.35 ± 0.02 respectively. 
Triangles are shifted by 0.1 T for clarity. 
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be the transition temperature Tj. For r < Ti, Fa should 
jump to zero which is not completely reproduced. We 
note that since rj G is not well defined in the isotropic fluid, 
it becomes problematic to extract Fa from eqn. Q be- 
low Ti. At r m , at the hexatic — * crystalline transition, 
Fa must diverge which indeed it does. This divergence 
can be identified with the divergence of the square of the 
translational correlation length £ + 0] 



fU(r)/fc B T ~ £ ~ exp (-— 



2c 



|i/r-i/r m | 



(•5) 



where c is again a constant and v — 0.36963. Fitting 
the values of Fa to the expression in eqn. © in the 
range 57.5 < T < 61 we obtain v = 0.35 ± 0.02 and 
r m = 61.3 ± 0.4 as an upper threshold. Again triangles 
represent evaluation of our data in sub-windows of 
variable size (same sizes as above). The finite size effect 
for r = 57.0 is negligible. Close to T m it increases but 
the values saturates for T — 59.1 and T = 60.8 and 
remain within the error-bars for the biggest sub- windows. 

In conclusion, we have checked quantitatively the 
change of quasi-long-range to short-range orientational 
order and extracted the correlation length i n the 
isotropic fluid and Frank's constant Fa in the hexatic 
phase from trajectories of a 2d colloidal system. We find 
a hexatic <-> isotropic liquid transition at Ti = 57.5 ±0.5. 
Three observations support this result: (i) the change of 
the distance dependence of Ge(r) (Fig. J3J), (ii) the con- 
dition FA(Ti) — 72/-7T for Frank's constant and (hi) the 
divergence of ^q. F° r the transition hexatic <-> crystal 
Fa diverges at T m . Both divergencies (extracted from 
just one correlation function) lead to critical exponents 
that are in good agreement with the KTHNY-theory. 
The measurements for melting and freezing support each 
other; so we may conclude that there is no hysteresis ef- 
fect of the phase-transitions. At the two transitions, the 
order parameters are observed to change continuously 
(within the resolution of T cx 1/T); no indication of a 
phase-separation (as for example stron g fl uctuations of 
the order parameters) has been found [3(j as has been 
reported by [17, 21]. So we believe that in our system - 
having a well-defined, purely repulsive pair-potential and 
a confinement to 2D that is free of any surface roughness 
- the transitions are second order. 

In [3ll l32| we verified that the Young's modulus 
becomes 167r at T m . We have now checked that Fa 
takes the value 72/-7T at Ti. These two findings together 
confirm the two-stage KTHNY melting scenario with its 
underlying microscopic picture of breaking the transla- 
tional symmetry by dislocation-pair- and orientational 
symmetry by disclination-pair-unbinding. 
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